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Abstract. The physics of gluon saturation and non-linear evolution at small values of parton momentum fraction x 
in the proton and nucleus is discussed in the context of experimental results at HERA and RHIC. The rich physics 
potential of low-x QCD studies at the LHC is discussed and some measurements in pp, pA and AA collisions accessible 
with the Compact Muon Solenoid (CMS) experiment are presented. 

PACS. 12.38.-t-24.85.+p-25.75.-q 



1 Introduction 

1.1 Parton structure and evolution 

Deep inelastic scattering (DIS) electron-proton ep (and 
electron-nucleus, eA) collisions provide a precise means to 
study the partonic structure of the proton (and nucleus). The 
inclusive DIS hadron cross section, d 2 a/dxdQ 2 , is a function 
of the virtuality Q 2 of the exchanged gauge boson (i.e. its "re- 
solving power"), and the Bjorken-jc fraction of the total nucleon 
momentum carried by the struck parton. The differential cross 
section for the neutral-current (y, Z exchange) process is writ- 
ten in terms of the target structure functions as 



d 2 a 
dxdQ 
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where Y± = 1 ± (1 — y) 2 is related to the inelasticity ij of the 
collision, and the structure functions F2,3,l(x,Q 2 ) describe the 
density of quarks and gluons in the hadron 1 . F 2 is the dom- 
inant contribution to the cross section over most of phase 
space. One of the most significant discoveries at HERA is the 
strong growth of the inclusive DIS cross section for decreas- 
ing Bjorken-x at fixed Q 2 as well as for increasing Q 2 at fixed 
x (Fig. 0. The strong scaling violations evident at small x in 
Fig.^are indicative of the increasing gluon radiation. At small 
x, F2 is sensitive to the sea quark distribution, driven by the 
gluon splitting, and since 



dF 2 (x/2,Q 2 )/dlnQ 2 - a s (Q 2 )xg(x,Q 2 ) , 



(2) 



the gluon density xg(x,Q 2 ) can be thus determined (Fig. |2j. 
Once measured at an input scale Qq > 2 GeV 2 , the parton 
distribution functions (PDFs) at any other Q 2 are given by 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evo- 
lution equations which govern the probability of parton branch- 
ings (gluon splitting, q,g— strahlung) in QCD |2|. 

1 F 2 °= e 1 x Li(qi + qi), xFt, °c xZ, (g,' -qi),F L <x U s xg. 
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Fig. 1. F 2 {x,Q 1 ) measured in DIS at HERA and fixed-target experi- 
ments. 



The DGLAP parton evolution, however, only takes into ac- 
count the Q 2 -dependence of the PDFs, resumming over single 
logarithms in a s \n{Q 2 ), "leading twist", but neglecting the 1 jx 
terms. At large energies (small x), the probability of emitting 
an extra gluon increases as °c a s ]n(l/x). In this regime, the 
evolution of parton densities proceeds over a large rapidity re- 
gion, Ay ~ ln(l/x), and the finite transverse momenta of the 
partons become increasingly important. Thus, their appropri- 
ate description is in terms of kj-unintegrated PDFs, xg(x,kj), 
described by the Balitski-Fadin-Kuraev-Lipatov (BFKL) equa- 
tion which governs parton evolution in x at fixed Q 2 [ 3 1 . 
Hints of extra BFKL radiation have been recently found in 
the enhanced production of forward jets at HERA compared 
to DGLAP expectations 1415 1 . 



David d'Enterria for the CMS collaboration: Low-x QCD with CMS at the LHC 
Hl+ZEUS Y = ln 1/x 



M 



20 



IS 



10 



Q = 20 GeV 2 ZEUS-JETS Fit 

| total uncert. 
HI PDF 2000 
Q 2 = 200GeV 2 ^exp. uncert. 

I total uncert. 




Fig. 2. Gluon distributions extracted at HERA (HI and ZEUS) as a 
function of x in three bins of Q 2 Q. 



1.2 Parton saturation and non-linear evolution at Iow a 

As shown in Figures ^ and |2] the gluon density rises very fast 
for decreasing x. For x < 0.01, the growth in F2 is well de- 
scribed by F 2 (x,Q 2 ) oc ^-^(e 2 ) with A, w 0.1 - 0.3 logarith- 
mically rising with Q 2 |6|. Eventually, at some small enough 
value of x one expects to enter a regime where the gluon density 
becomes so large that non-linear (gg fusion) effects become 
important, taming the growth of the parton densities. In such a 
high-gluon density regime three things are expected to occur: 
(i) the standard DGLAP and BFKL linear equations should no 
longer be applicable since they only account for single par- 
ton branchings (1 — >2 processes) but not for non-linear (2 — >1) 
gluon recombination; (ii) pQCD (collinear and kj) factoriza- 
tion should break due to its (now invalid) assumption of inco- 
herent parton scattering; and, as a result, (iii) standard pQCD 
calculations lead to a violation of unitarity even for Q 2 > A 2 . 
Figure[5]schematically depicts the different domains of the par- 
ton density as a function of y = ln(l jx) and Q 2 . The transition 
to the regime of saturated PDFs is expected for small x val- 
ues below an energy-dependent "saturation momentum", Q s , 
intrinsic to the {size of the) hadron. Since xg(x, Q 2 ) can be inter- 
preted as the number of gluons with transverse area r 2 ~ 1 /Q 2 
in the hadron wavefunction, an increase of Q 2 effectively di- 
minishes the 'size' of each parton, partially compensating for 
the growth in their number (i.e. the higher Q 2 is, the smaller 
the x at which saturation sets in). Saturation effects are, thus, 
expected to occur when the size occupied by the partons be- 
comes similar to the size of the hadron, TiR 2 . This provides a 
definition for the saturation scale of an arbitrary hadron with A 
nucleons (i.e. with gluon density xG = A ■ xg): 



1 , 



In CT(Y)=AY 




In Q' 



Fig. 3. QCD "phase diagram" in the 1 fx, Q 2 plane (each dot repre- 
sents a parton with transverse area ~ 1 /Q 2 carrying a fraction x of the 
hadron momentum) (7J. 

with A. sa 0.25 1 8 1. Equation Q indicates that Q s grows with the 
number of nucleons, A, of the target, and the energy of the colli- 
sion, yfs, or equivalently, the rapidity of the gluon y — ln( 1 jx) . 
The mass number, A, dependence implies that, at equivalent en- 
ergies, saturation effects will be enhanced by factors as large as 
A 1 / 3 w 6 in heavy nuclear targets (A = 208 for Pb) compared 
to protons. In the last fifteen years, an effective field theory of 
QCD in the high-energy (high density, small x) limit has been 
developed - the Colour Glass Condensate (CGC) J5) - which 
describes the hadrons in terms of classical fields (saturated 
gluon wavefunctions) below the saturation scale Q s . The sat- 
uration momentum Q s introduces a (semi-)hard scale, Q s A, 
which not only serves as an infrared cut-off to unitarize the 
cross sections but allows weak-coupling perturbative calcula- 
tions (0Cj(<2i) -Cl) in a strong colour field background. 
Hadronic and nuclear collisions are seen as collisions of classi- 
cal wavefunctions which "resum" all gluon recombinations and 
multiple scatterings. The quantum evolution in the CGC ap- 
proach is given by the JIMWLK 1 10 1 non-linear equations (or 
by their mean-field limit for N c — > °°, the Balitsky-Kovchegov 
equation 1 1 1 1) which reduce to the standard BFKL kernel at 
higher x values. 



2 Parton saturation: Experimental studies 

The main source of information on the quark densities is ob- 
tained from measurements of (i) the structure functions F% 3 in 
lepton-hadron scattering, and (ii) lepton pair (Drell-Yan) pro- 
duction in hadron-hadron collisions. The gluon densities, xG 
enter at LO directly in hadron-hadron scattering processes with 
(i) prompt photons and (ii) jets in the final state, as well as in 
the (difficult) measurement of (iii) the longitudinal DIS struc- 
ture function Fi (and also indirectly in F2 through the derivative 
in Eq. 0). In addition, (iv) heavy vector mesons (7/\|/, Y) from 
diffractive photoproduction processes 2 are a valuable probe of 
the gluon density since their cross sections are proportional to 



,1/3 



(3) 



2 Diffractive jp (yA) processes are characterized by a quasi-elastic 
interaction - mediated by a Pomeron or two gluons in a colour singlet 
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Fig. 4. Experimental measurements at various facilities providing in- 
formation on the gluon PDF in different ranges of Bjorken-x. 



non-perturbative range (Q 2 < 1 GeV 2 ) dominated by Regge 
dynamics rather than quark/gluon degrees of freedom. An ex- 
ample of the current lack of knowledge of the nuclear densi- 
ties at low x is presented in Fig. [6] where different available 
parametrizations of the ratio of Pb to proton gluon distribu- 
tions, consistent with the available nDIS data at higher x, show 
differences as large as a factor of three 11191 . 




the square of xG 1121131 : 



da 



ypA^VpA 



dt 



asYee \^[xG{x,Q 2 )]\ (4) 
with Q 2 = M 2 >/4 and x=M^/W 2 pA . (5) 



(=0 



3aM 5 v 



In hadronic collisions, one commonly measures (real and 
virtual) photons and jets at central rapidities (y = 0) where 
x — xj — Q/ y/s, with Q ~ pr,M the characteristic scale of 
the hard scattering. However, one can probe smaller xi values 
in the target by measuring the corresponding cross sections in 
the forward direction. Indeed, for a 2 — * 2 parton scattering the 
minimum momentum fraction probed in a process with a parti- 
cle of momentum pr produced at pseudo-rapidity T| is 1 14 1 



min *T € 

X 2 



-i 



2-x T e^ 



where xj = 2pr / \fs , 



(6) 
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Fig. 5. Available measurements in the (x, Q 2 ) plane used for the deter- 
mination of the proton 1 16 1 (top) and nuclear 1 17 1 (bottom) PDFs. 



i.e. x%'" decreases by a factor of ^10 every 2 units of rapidity. 
Though Eq. is a lower limit at the end of phase-space (in 
practice the (x2) values in parton-parton scatterings are at least 
10 larger than x'™ 11 1141 1. it provides the right estimate of the 
typical X2 — (pr / V^) e ^ values reached in non-linear 2^1 
processes (in which the momentum is balanced by the gluon 
"medium") as described in parton saturation models 1151 . 

Figure 0] summarizes the range of experimental processes 
sensitive to the gluon density and their approximate x cover- 
age. Figure|5]shows the kinematical map in (x,Q 2 ) of the DIS, 
DY, direct y and jet data used in the PDF fits. Results from 
HERA and the Tevatron cover a substantial range of the proton 
structure (1(T 4 < x < 0.8, 1 < Q 2 < 10 5 GeV 2 ) but the avail- 
able measurements are much rarer in the case of nuclear tar- 
gets (basically limited to fixed-target studies, 10~ 2 < x < 0.8 
and 1 < Q 2 < 10 2 GeV 2 ). As a matter of fact, the nuclear parton 
distributions are basically unknown at low x (x < 0.01) where 
the only available measurements are fixed-target data in the 



state - in which the p (A) remains intact (or in a low excited state) and 
separated by a rapidity gap from the rest of final-state particles. 



2.1 HERA results 

Though the large majority of ep DIS data collected during 
the HERA-I phase are consistent with standard DGLAP pre- 
dictions, more detailed and advanced experimental and the- 
oretical results in the recent years have pointed to interest- 
ing hints of non-linear QCD effects in the data. Arguably, 
the strongest manifestation of such effects is given by the so- 
called "geometric scaling" property observed in inclusive Odis 
for x < 0.01 1 20 1 as well as in various diffractive cross sec- 
tions 12 1 1221 . For inclusive DIS events, this feature manifests 
itself in a total cross section at small x (x < 0.01) which is 
only a function of x = Q 2 / Q 2 {x), instead of being a function 
of x and Q 2 /Q 2 separately (Fig. 0. The saturation momen- 
tum follows Q s (x) = Qo(x/xq)^ with parameters A, ~ 0.3, Qo 
= 1 GeV, and xq ~ 3T0~ 4 . Interestingly, the scaling is valid 
up to very large values of x, well above the saturation scale, in 
an "extended scaling" region with Q 2 < Q 2 < 2 4 /Ag CD 1231 
0. The saturation formulation is suitable to describe not only 
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Fig. 6. Ratios of the Pb over proton gluon PDFs versus x from different 
models at Q 2 = 5 GeV 2 . Figure taken from 1181 . 
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Fig. 7. Geometric scaling in the DIS J*p cross sections plotted versus 
X = Q 2 /Q 2 in the range x < 0.01, 0.045 < Q 2 < 450 GeV 2 l20l . 



inclusive DIS, but also inclusive diffraction y*p — > X p. The 
very similar energy dependence of the inclusive diffractive and 
the total cross section in y*p collisions at a given Q 2 is eas- 
ily explained in the Golec-Biernat Wusthoff model |20| but 
not in standard collinear factorization. Furthermore, geomet- 
ric scaling has been also found in different diffractive DIS 
cross sections (inclusive, vector mesons, deeply-virtual Comp- 
ton scattering DVCS) 1 12 11221 . All these results suggest that 
the observed scalings are indeed manifestations of the satura- 
tion regime of QCD. Unfortunately, the value of Q s ~ 1 GeV 
at HERA lies in the transition region between the soft and hard 
sectors and, therefore, non-perturbative effects obscure the ob- 
tention of clearcut experimental signatures. 



2.2 RHIC results 

The expectation, based on Eq. 0, of enhanced parton satura- 
tion effects in the nuclear wave functions accelerated at ultra- 
relativistic energies has been one of the primary physics moti- 
vations for the heavy-ion programme at RHIC 3 1 9 1 . Further, the 
properties of the high-density matter produced in the final-state 
of AA interactions cannot be properly interpreted without hav- 
ing determined the influence of initial state modifications of 
the nuclear PDFs. In this context, after five years of operation, 
two main experimental observations at RHIC have been found 
consistent with CGC predictions: (i) the modest hadron multi- 
plicities measured in AuAu reactions, and (ii) the suppressed 
hadron yield at forward rapidities in dAu collisions. 

The bulk, dN c h/dr[\^=o ~ 700, multiplicities measured at 
mid-rapidity in central AuAu at ■Js^, = 200 GeV are compar- 
atively lower than the dN c h / dr\\^=o ~ 1000 predictions |24| of 
"minijet" scenarios, soft Regge models, or extrapolations from 
an incoherent sum of proton-proton collisions, but they can 
be reproduced by approaches based on gluon saturation 181251 
which take into account a reduced parton flux in the nuclear tar- 
gets, i.e. f a /A(x, Q 2 ) <A- f a /N(x,Q 2 ). In the CGC calculations, 
the final hadron multiplicities are assumed to be simply related 
to the initial number of released partons (local parton-hadron 
duality) which are depleted in the initial state compared to pp 
collisions due to non-linear gluon-gluon recombinations |8|. 
Simple assumptions related to the dependence of the saturation 
scale on energy and overlapping area of the colliding nuclei, 
describe the centrality and center-of-mass (cm.) energy depen- 
dences of the bulk AA hadron production (Fig. [8}. 

Vs (GeV/A) 




Fig. 8. Dependences on cm. energy and centrality (given in terms 
of the number of nucleons participating in the collision , iVp art) of 
dN c h/dl]\^ = o (normalized by Np^): PHOBOS AuAu data |26| vs the 
predictions of the saturation approach 1251 . 

The second manifestation of saturation-like effects in the 
RHIC data is the BRAHMS observation |^) of suppressed 
yields of moderately high-pr hadrons (pr ~2-4 GeV/c) in 
dAu relative to pp collisions at r| w 3.2. Hadron production at 

3 The saturation scale at y =0 in Au at RHIC is Q 2 ~ 2 GeV, much 
larger than that of protons probed at HERA, Q 2 ~ 0.5 GeV. 
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Fig. 9. Nuclear modification factor R^au (Pt ) f° r charged hadrons pro- 
duced in dAu at = 200 GeV: BRAHMS data [ 27 1 versus DGLAP 
shadowing 1 15 1 and CGC 1 28 1 predictions. Fig. adapted from 1151 . 

such small angles is sensitive to partons in the Au nucleus with 
x w i^(10~ 3 ). The observed nuclear modification factor, RdAu ~ 
0.8, cannot be reproduced by pQCD calculations that include 
standard leading-twist shadowing of the nuclear PDFs 11511 41 
but can be described by CGC approaches 1 28 1 that parametrize 
the Au nucleus as a saturated gluon wavefunction. In addition, 
a recent analysis of the nuclear DIS Fi data also confirms the 
existence of "geometrical scaling" for* <0.017 |25|. 



3 Low-jc QCD at the LHC 




x 

Fig. 10. Kinematical (x, Q 2 ) range probed at various rapidities y and 
cm. energies in = 8.8 TeV pA collisions at the LHC 1181 . 

(iii) Complete electromagnetic (EM) and hadronic (HAD) 
calorimetry for full jet reconstruction over |t|| < 3 and A0 
= 2k with a large statistical significance for single jet and 
jet+X (X = jet, y, Z) channels. 

(iv) Unparalleled forward physics capabilities thanks to the 
forward hadronic calorimeter (HF, 3< |t|| <5), TOTEM 
Tl (3.1 < |r|| < 4.7) and T2 (5.5< |r)| <6.6) trackers, and 
CASTOR (5.3 < |r|| <6.7) and zero-degree (ZDC, |r)| >8.1 
for neutrals) calorimeters. 



The Large Hadron Collider (LHC) at CERN will provide pp, 
pA and AA collisions at y/s^ = 14, 8.8 and 5.5 TeV respec- 
tively with luminosities «5f ~ 10 34 , 10 29 and 5 TO 26 cm -2 
s _1 . Such large cm. energies and luminosities will allow de- 
tailed QCD studies at unprecedented low x values thanks to 
the copious production of hard probes (jets, quarkonia, prompt 
y, Drell-Yan pairs, etc.). The expected advance in the study of 
low-x QCD phenomena will be specially substantial for nuclear 
systems since the saturation momentum, Eq. (0, Q 2 . « 5 - 10 
GeV 2 , will be in the perturbative range |8 1, and the relevant x 
values, Eq. 0, will be 30-70 times lower than AA and pA re- 
actions at RHIC: x » 10 3 (10 5 ) at central (forward) rapidities 
for processes with Q 2 -10 GeV 2 (Fig.fTol. 



3.1 The CMS experiment 

The CMS experiment is one of the two large general-purpose 
detectors being installed at the LHC. Its experimental capabili- 
ties are extremely well adapted for the study of low-x phenom- 
ena with proton and ion beams featuring: 

(i) Very large acceptance at midrapidity (|t|| < 2.5, full (])) for 
charged and neutral hadrons as well as fi , e , and y over a 
wide range of pj (the 4 T magnetic field results in the best 
track momentum resolution at LHC). 

(ii) Excellent muon reconstruction leading to the best mass 
resolution for J/y, and T measurements at the LHC. 




Fig. 11. Layout of the detectors in the CMS forward region. 

The combination of HF, TOTEM, CASTOR and ZDC 
(Fig. Hit makes of CMS the largest acceptance detector ever 
built at a hadron collider. The HF |29|, located 11.2 m away 
on both sides of the interaction point (IP), is a steel plus 
quartz-fiber Cerenkov calorimeter with 1200 channels (Ar| x 
A(j) —0.18x0.18, 1.65 m absorber corresponding to 10. 3A./) 
sensitive to the deposited EM and HAD energy, allowing 
jet reconstruction at very forward rapidities. The Tl and 
T2 telescopes are part of the the TOTEM experiment 1 30 1 
which shares the same IP as CMS and are mainly designed 
to measure charged tracks from diffractive dissociation pro- 
cesses. CASTOR [31 1 is an azimuthally symmetric electromag- 
netic/hadronic calorimeter situated at 14.37 m from the interac- 
tion point covering the same acceptance as T2. The calorime- 
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ter is a Cerenkov-light device, consisting of successive layers 
of tungsten absorber and fused silica (quartz) plates as active 
medium arranged in 2 EM (10 HAD) sections of about 22Xq 
(10. 3X]) radiation (interaction) lengths. The ZDC 1321 is also 
a tungsten+quartz sampling Cerenkov calorimeter with 5 EM 
(19Xq, divided in x) and 4 HAD (5.6A-/, divided z) sections. It is 
located at 140 m from the CMS vertex at the end of the straight 
sections of the two LHC pipes containing the countercirculat- 
ing beams. The purpose of the ZDC is to measure very forward 
going neutrons and photons with ~10% (2 mm) energy (posi- 
tion) resolution. 



3.2 Low-x QCD measurements in CMS 

The following three measurements in pp, pA and AA collisions 
are being considered in CMS to look for signatures of high 
gluon density effects at low x. 



3.2.2 QQ photoproduction (electromagnetic AA collisions) 

High-energy diffractive production of heavy vector mesons 
(7/\|/,T) proceeds through colourless two-gluon exchange 
(which couples to y — > QQ) and is thus a sensitive probe of the 
low-x gluon densities, see Eq. @. Ultra-peripheral interactions 
(UPCs) of high-energy heavy ions generate strong electromag- 
netic fields which help constrain the low-x behaviour of xG via 
quarkonia produced in y-nucleus collisions 1 34 1 . Lead beams at 
2.75 TeV have Lorentz factors y = 2930 leading to maximum 
(equivalent) photon energies co mflA «y/R~ 100 GeV, and cm. 
energies W™' x w 160 GeV and W$' x » 1 TeV. From Eq. ©, 
the x values probed in yA — * J/\\tA processes at y = 2 can be as 
low as x <~ 10~ 5 . The CMS experiment can measure T — > 
produced in electromagnetic PbPb collisions tagged with neu- 
trons detected in the ZDCs (as done at RHIC 1341 ). Figure 13 
shows the expected dimuon invariant mass distributions pre- 
dicted by STARLIGHT |35) within the CMS acceptance for an 
integrated PbPb luminosity of 0.5 nb~ 1 . An T peak with ~ 1200 
counts is clearly seen on top of the continuum. 



3.2.1 Forward jets (pp, pA, AA) 



The cross section for dijet production in the forward direc- 
tion, "Miiller-Navelet" jets |33|, is a particularly sensitive mea- 
sure of the small x parton dynamics in hadronic collisions |5 1. 
The two HF calorimeters (3< |r|| <5), specifically designed to 
measure energetic forward jets 4 , have an energy resolution of 
—20% for typical jets with E T — 40 GeV (i.e. E=E T coshT) « 
1 TeV at T|= 4). In the presence of low-x saturation effects, the 
forward-backward dijet production cross section (separated by 
Ar| — 9 and, thus, measurable in each of the HFs) is expected 
to be suppressed by a factor of —2 in pp at 14 TeV (Fig. I12> . 
A study is underway to determine the feasibility of such mea- 
surements in CMS. 




Q (GeV) 



Fig. 12. Ratio of the saturation over BFKL predictions for the Miiller- 
Navelet forward jet cross sections in pp collisions at y/s = 14 TeV as 
a function of Q = Q\ = Q2 for different values of y = y\ = — yi |5|. 
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Fig. 13 
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invariant mass from yPb — > YPb* — > 
as given by Starlight 03 for UPC PbPb 
5.5 TeV in the CMS acceptance. 



■H + H 



3.2.3 Forward Drell-Yan pairs (pp, pA, AA) 

High-mass Drell-Yan pair production at the very forward ra- 
pidities covered by CASTOR and T2 (|t|| — 5 — 6) can probe 
the parton densities down to x — 10~ 6 . A study is underway in 
CMS to combine the CASTOR electromagnetic energy mea- 
surement together with the good position resolution of T2 for 
charged tracks, to trigger on and reconstruct the e + e~ invariant 
mass in pp collisions at 14 TeV, and perform a two-dimensional 
study of xg in the M 2 and x plane. 



4 Conclusion 



4 The HF plays a prominent role in forward jet tagging for the 
vector-boson-fusion (qq — > qqH) Higgs production channel. 



We have reviewed the physics of non-linear QCD and high 
gluon densities at small fractional momenta x with emphasis 
on the existing data at HERA (proton) and RHIC (nucleus) 
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which support the existence of a parton saturation regime (also 
known as colour-glass-condensate). The future perspectives at 
the LHC have been presented, including the promising capabil- 
ities of the forward CMS detectors to study the parton densities 
down to x r-j 10~ 6 with various hard probes (jets, quarkonia, 
Drell-Yan). The programme of investigating the dynamics of 
low-x QCD is not only appealing in its own right but it is an 
essential prerequisite for predicting a large variety of hadron- 
, photon- and neutrino- scattering cross sections at very high 
energies. 



Acknowledgments 

The author acknowledges valuable comments from R. Vogt and 
J. Jalilian-Marian, and thanks J. Nystrand for providing the 
STARLIGHT photoproduction cross sections in PbPb at LHC. 
This work is supported by the 6th EU Framework Programme 
contract MEIF-CT-2005-025073. 



24. K. J. Eskola, Nucl. Phys. A 698 (2002) 78. 

25. N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. 
Lett. 94 (2005) 022002. 

26. B. B. Back et al. [PHOBOS], Nucl. Phys. A757 (2005) 28. 

27. I. Arsene et al. [BRAHMS], Nucl. Phys. A757 (2005) 1. 

28. D.Kharzeev, Y.Kovchegov, K.Tuchin, Phys. Lett. B599 (2004) 23 

29. A. S. Ayan et al, J. Phys. G 30 (2004) N33. 

30. V. Berardi et al. [TOTEM], "TOTEM: Technical design report", 
CERN-LHCC-2004-002. 

31. P. Katsas et al, submitted to Nucl. Instrum. Methods A. 

32. O. A. Grachov et al., Proceeds. CALOR'06. lrmcT- ex/0608052 

33. A. H. Mueller and H. Navelet, Nucl. Phys. B 282 (1987) 727. 

34. D. d'Enterria, Proceeds. Quark Matter' 05, nucl-ex/0601001 

35. S. R. Klein, J. Nystrand, Phys. Rev. C60 (1999) 014903; A. Baltz, 
S. Klein, J. Nystrand, Phys. Rev. Lett. 89 (2002) 012301. 



References 

1. M. Dittmar et al, in Proceeds. HERA and the LHC, 
hep-ph/0511119 

2. V.N. Gribov and L.N. Lipatov. Sov. Journ. Nucl. Phys. 15 (1972) 
438; G. Altarelli and G. Parisi, Nucl. Phys.B126 (1977) 298; Yu. 
L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

3. L.N. Lipatov, Sov. J. Nucl. Phys. 23 (1976) 338; E.A. Kuraev, 
L.N. Lipatov and V.S. Fadin, Zh. Eksp. Teor. Fiz 72, (1977) 3; 
Ya.Ya. Balitsky, L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

4. S. Chekanov et al. [ZEUS Collab.], Phys. Lett. B 632 (2006) 13; 
A. Aktas et al. [HI Collab.], Eur. Phys. J. C 46 (2006) 27. 

5. C. Marquet and C. Royon, Nucl. Phys. B 739 (2006) 131. 

6. C. Adloff et al. [HI Collab.], Phys. Lett. B 520 (2001) 183. 

7. E. Iancu, hep-ph/0608086 

8. D. Kharzeev and M. Nardi, Phys. Lett. B 507 (2001) 121; 
D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. A 747 (2005) 609. 

9. See e.g. E. Iancu and R. Venugopalan, in QGP. Vol 3 Eds: R.C. 
Hwa and X.N. Wang, World Scientific, Singapore, hep-ph/0303204 
J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. Nucl. Phys. 56 
(2006) 104; andrefs. therein. 

10. J. Jalilian-Marian, A. Kovner, A. Leonidov and H. Weigert, Nucl. 
Phys. B504 (1997) 415; Phys. Rev. D59 (1999) 014014; E. Iancu, 
A. Leonidov and L. McLerran, Nucl. Phys. A692 (2001) 583; 

1 1. I. Balitsky, Nucl. Phys. B 463, (1996) 99; Yu.V. Kovchegov, Phys. 
Rev. D 61, (2000) 074018. 

12. M. G. Ryskin et al., Z. Phys. C76, (1997) 231. 

13. T. Teubner, Proceeds. DIS'05, AIP Conf. Proc. 792 (2006) 416. 

14. V.Guzey, M.Strikman, W.Vogelsang, Phys. Lett. B603 (2004) 173 

15. A. Accardi, Acta Phys.Hung. A22 (2005) 289. 

16. P. Newman, Int. J. Mod. Phys. A 19 (2004) 1061. 

17. D. d'Enterria, J. Phys. G 30 (2004) S767. 

18. A. Accardi et al, in CERN Yellow report on Hard probes in 
Heavy Ion collisions at the LHC, hep-ph/0308248 

19. N. Armesto, J. Phys. G to appear, hep-ph/0604108 

20. K. Golec-Biernat and M. Wusthoff, Phys. Rev. D59, (1999) 
014017; Phys. Rev. D60, (1999) 114023. 

21. J. R. Forshaw and G. Shaw, JHEP 0412 (2004) 052. 

22. C. Marquet, L. Schoeffel, hep-ph/0606079 

23. E. Iancu, K. Itakura, L. McLerran, Nucl. Phys. A708 (2002) 327. 



